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The 7 -ray binary LS I +61°303 shows a discontinuity of the periodicity in its GeV emission. In this paper, we 
show that during the epochs when the timing analysis fails to determine the orbital periodicity, the periodicity 
is in fact present in the two orbital phase intervals ^ = 0.0 — 0.5 and ^ = 0.5 — 1.0. That is, there are two 
periodic signals, one towards periastron (i.e., <1> = 0.0 — 0.5) and another one towards apastron (<1> = 0.5 — 1.0). 
The apastron peak shows the same orbital shift as the radio outburst and, in addition, reveals the same two 
periods Pi and P 2 that are present in the radio data. The 7 -ray emission of the apastron peak normally just 
broadens the emission of the peak around periastron. Only when it appears at ^ = 0.8 — 1.0 because of the 
orbital shift, it is detached enough from the first peak to become recognizable as a second orbital peak, which 
is the reason why the timing analysis fails. Two 7 -ray peaks along the orbit are predicted by the two-peak 
accretion model for an eccentric orbit that was proposed by several authors for LS I +61° 303. 


1. Introduction 

The stellar system LS I +61°303 is a member of 
the small class of 7 -ray binaries, which are defined 
as binary stars with a peak in the spectral energy 
distribution above 1 MeV [T]. A sketch of the system 
is shown in Fig. LS I +61°303 consists of a Be star 
and a compact object in an eccentric orbit, e = 0.72± 
0.15 [5], with orbital period Pi = 26.4960 ± 0.0028 d 
[ 6 ]. The Be star is rapidly rotating and losing mass 
in form of an equatorial disk [5]. The nature of the 
compact object could not yet be established, because 
the masses are poorly constrained due to the large 
uncertainty in the inclination angle [5]. The orbital 
phase of the binary system is defined as 

*=+ - (+) • 

where to = MJD 43366.275 [ 6 ]. Periastron occurs at 
orbital phase 4> = 0.23 [5]. 

Radio outbursts are observed at orbital phases 4> = 
0.5 — 0.9, i.e., around apastron. Their peak flux den¬ 
sities are modulated in amplitude and orbital phase 
occurrence by a long-term period Piong = 1667 ± 8 d 
[ 6 ]. The long-term phase 0 is defined analogous to 
the orbital phase 4> by replacing Pi by Piong- 

The source LS I +61°303 is highly variable and peri¬ 
odic all over the electromagnetic spectrum from radio 
to very high energy 7 -rays [U [SJ El [7] . The GeV 7 -ray 
light curve, as obtained using Fermi LAT data, has so 
far been reported to peak at orbital phases around pe¬ 
riastron mm (see Fig.[^. Timing analysis shows that 
the orbital period is present in the Fermi LAT light 
curve from this source, however not with equal power 
all of the time la nn]. There are times (0-phases) 
when the period is outstanding and there are times 
when the period is completely absent from the power 
spectrum, as shown well in Fig. 4 of HO]. Moreover, 
Fig. 3 of [ 10 ] shows that GeV data also show the long¬ 
term periodical variation affecting the radio data, but 
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Figure 1: Sketch of LS I +61° 303. The periodic GeV 
outburst at periastron was first reported by [2]- In 
Sect. 3 a newly discovered periodic apastron GeV peak is 
discussed [3]. The radio outburst occurs only at 
apastron, whereas at periastron, only a low level of 
emission is present [see Fig. 1-Right in [4]. 


only at a specific orbital phase interval, 4> = 0.5 — 1.0, 
that is around apoastron. 

We are aimed here to investigate the discontinuity 
in the periodicity of the GeV 7 -ray emission at peri¬ 
astron, the possible relationship of its disappearance 
with the variation of the emission around apastron, 
and finally the possible relationship between GeV and 
radio emission. 


2 . Data analysis 

For the present analysis |3] we use Fermi LAT 
data from LS I +61°303 spanning the time range 
August 5, 2008 (MJD 54683) until June 30, 2014 
with an energy range of 100 MeV to 300 GeV. 
For the computation of the light curves we used 
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the script like_lc.pl written by Robin Corbet. ^ 
Only source-event-class photons were selected for 
the analysis. Photons with a zenith angle greater 
than 100° were excluded to reduce contamination 
from the Earth’s limb. For the diffuse emission we 
used the model glLiem_v05_revl.fit and the template 
iso_source_v05_revl.txt. We used the instrument re¬ 
sponse function (IRE) P7REP/background_revl, and 
the model file was generated from the 2FGL catalogue 
[TT], all sources within 10° of LS I H-61°303 were in¬ 
cluded in the model. LS I -1-61°303 was fitted with a 
log-parabola spectral shape and with all parameters 
left free for the fit, performing an unbinned maximum 
likelihood analysis. The other sources were fixed to 
their catalogue values. We produced light curves with 
a time bin size of one day and of five days. 

We investigated [for details see |3] the temporal 
evolution of the orbital periodicity by means of a 
wavelet analysis [12] and Lomb-Scargle timing anal¬ 
ysis mM- 

3. Results: A periodic signai around 
apastron 

3.1. Waveiet anaiysis 

Our results are shown Fig. 2. The first plot of Fig. 2 
presents the examined data set. The wavelet analy¬ 
sis was applied to the y-ray data vs time, however, 
for a straightforward comparison with radio data, we 
express in the other plots of Fig. 2 the x-axis as the 
long-term phase 0. The second plot of Fig. shows 
the wavelet plot for the whole data set, i.e., the whole 
orbital period $ = 0.0—1.0. The absence of the orbital 
period around 0 « 7.2 is consistent with the previous 
finding shown in Fig. 4 of m- When wavelet analy¬ 
sis is performed only on data from the orbital phase 
intervals 4> = 0.0 — 0.5 (middle) and <1> = 0.5 — 1.0 
(bottom), it is revealed that there is always a periodic 
signal at ^ = 0.0 —0.5 (periastron). Moreover, there is 
a periodic signal at ^ = 0.5 — 1.0 (apastron). The lat¬ 
ter becomes particularly strong during the time when 
the orbital period is absent from the power spectra of 
^ = 0.0-1.0 [3]. 

3.2. Lomb-Scargle timing analysis 

Figure shows Lomb-Scargle periodograms of the 
y-ray flux from LS I +61°303. The data have been 
selected from orbit phase intervals like in the previ¬ 
ous section. In the periodogram for the entire orbit 
(Fig. I a) the strongest feature is a peak which agrees 


^ http://f ermi.gsf c.nasa.gov/ssc/data/analysis/user/ 


well with the orbital period Pi found by [TT]. Fig- 
ures id, e, and f, refer to only data from 4> = 0.5 —1.0. 
In this orbital phase interval the peak at Piong is a very 
strong feature, in agreement with the findings of m- 
Moreover, the zoom of Fig. |3]d, i.e., Fig. [^e, shows a 
second peak, P 2 = 26.99 ± 0.08 d. This second peak 
becomes stronger and is almost as strong as the peak 
at Pi = 26.48 ± 0.08 d in the 5 day integrated data 
in Fig. [^f. The periods Pi, P 2 (see Fig. 4), and Piong 
here present are typical periodicities in radio data as 
shown in m- 

3 . 3 . Folded Fermi LAT data: The 
apastron GeV peak and its orbital shift 

Figure 5 a and b show Fermi LAT data from the 
time (0 ~ 7.2) of the disappearance of the orbital 
period from the power spectra folded with the orbital 
period. A second peak is evident at orbital phases 
4> = 0.8 — 1.0. Figure 5 c shows Fermi LAT data for 
another 0. It is now interesting to compare these plots 
with radio data. Because of the long-term periodicity 
we can compare y-ray and radio data having the same 
fractional part of 0. Figure 5d shows GBI radio data 
at 8 GHz [for details see|3]. 


4 . Conclusions 

During the intervals where the orbital periodicity is 
absent from the power spectra, wavelet and the folded 
light curves show two periodic signals, one at perias¬ 
tron and a second at apastron. The presence of the 
second periodic outburst disturbs the timing analysis 
and prevents it from finding the orbital periodicity. 
Gomparison with the folded radio data (Fig. 4d) sug¬ 
gests that the apastron GeV peak follows the same or¬ 
bital shift as the radio outbursts [3]. It is well-known 
the phenomenon of the orbital shift of the radio out¬ 
burst in LS I +61°303: The largest outbursts occur at 
orbital phase 0.6, afterwards, with the long-term pe¬ 
riodicity, the orbital phase of the peak of the outburst 
changes, as analysed by m in terms of orbital phase 
shift, by m in terms of timing residuals, and repro¬ 
duced recently by the precessing jet model in [18] , here 
shown in Fig. 5. 

Our result of two GeV peaks along the orbit corrob¬ 
orates the two-peak accretion model for LS I +61° 303. 
The hypothesis that a compact object that accretes 
material along an eccentric orbit undergoes two ac¬ 
cretion peaks along the orbit was suggested and devel¬ 
oped by several authors for the system LS I +61°303 
The first accretion peak is predicted 
to occur close to the Be star and to give rise to a ma¬ 
jor high-energy outburst. The second accretion peak 
is predicted to occur much farther away from the Be 
star, where the radio outburst occurs, and a minor 
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high-energy outburst is predicted there [21] . The pre¬ 
dicted periastron event corresponds well to the ob¬ 
served GeV peak towards periastron, the second pre¬ 
dicted high-energy outburst, corresponds well to the 
here discussed apoastron peak. 
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Figure 2: Wavelet analysis of Fermi-LAT data. The strength of periodicity is colour coded as indicated in the bottom 
bar. (a) Fermi-LAT data with a time bin of 1 d. (b) Wavelet analysis for the whole orbital interval 0.0 — 1.0 (b-d use a 
time bin of one day), (c) Wavelet analysis for the orbital interval 4> = 0.5 — 1.0, i.e., around apoastron. (d) Wavelet for 
eConfl]6)l^ft30ftJerval 4> = 0.0 — 0.5, i.e., around periastron. 
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Figure 3: Lomb-Scargle periodogram of the Fermi LAT data (with a time bin of one day). Figure 3 in [3]. (a) Data in 
the orbital phase $ = 0.0 — 1.0. (b) Zoom of Fig. |^a. (c) Same as |3]b for data with a time bin of 5 d. (d) Data in the 
orbital phase $ = 0.5 — 1.0. The periods P 2 and Piong here present are typical periodicities in radio data m- (e) Zoom 
of Fig. id. (f) Same as[^e for data with a time bin of 5 d. (g) Data in the orbital phase $ = 0.0 — 0.5. (h) Zoom of 
Fig. I^g. (i) Same as[^h for data with a time bin of 5 d. 
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Figure 4: (a)-(c) Folded Fermi LAT 7 -ray data (100 MeV - 300 GeV). The blue curve in (c) is that of [9]. (d) Folded 

GBI 8 GHz radio data. The here discovered periodic apastron GeV peak follows the same timing characteristic (i.e., Pi 
and P 2 are both present) as the periodic radio peak, which also occurs around apastron. During the time when the 
orbital period disappears from the power spectra of the 7 -ray light curve (see Fig. 2 b, © 7.2) the apastron GeV peak 

becomes well visible in the folded light curve, because it is more displaced from the periastron peak (see [3] and here 
Sect. 4). 



Figure 5: Left: Timing analysis of 6.7 years of GBI radio data at 2 and 8 GHz results in two periods. 

Pi = 26.49 ± 0.07d, P 2 = 26.92 ± 0.07d. The long-term period Piong = 1667 ± 8 d is consistent with the period 
Pbeat = l/(z^i — n 2 ) = 1667 ± 393 d resulting from the beating between the two close periodicities Pi and P 2 [E]. Right: 
Orbital shift of the radio outburst of LS I +61°303 in the precessing jet model of [18]. At © = 0.86 the outbursts peak 
at 4> 0.6. At © = 0.265 the outbursts peak at 4> 0.85. 
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